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PREFACE 


During  its  Meeting  in  Aix-en-Provence,  in  Fall  1980,  a  paper  was  presented  by 
Mr  H.Gfidel  to  the  Sub-Committee  on  Aeroelasticity  of  the  Structures  and  Materials  Panel 
on  the  use  of  structural  optimization  methods  to  obtain  practical  minimum-weight  designs 
which  meet  the  constraints  of  loads,  stresses,  buckling,  deflections,  divergence  and  flutter. 
Impressive  results  were  exhibited. 

The  excellent  presentation  of  this  advanced  application  of  optimization  procedures 
was  welcomed  by  the  Sub-Committee  as  being  a  very  important  contribution  to  its  activities 
in  Aeroelasticity.  Its  publication  as  an  AGARD  Report  will  help  dissemination  of  new 
methods  and  techniques  amongst  the  NATO  community. 


G.COUPRY 

Chairman,  Sub-Committee 
on  Aeroelasticity 


APPLICATION  OF  A  STRUCTURAL  OPTIMIZATION  PROCEDURE 
FOR  ADVANCED  WINGS 
by 

H.  G6del 

Aerospace  Engineer,  Aeroelastic  Branch 
and 

G.  Schneider 

Aerospace  Engineer,  Aeroelastic  Branch 

MESSERSCHMITT-BflLKCW-BLOHM  GMBH 
Airplane  Division 
P.O.Box  80  11  60,  8000  Munich  80 
W. -Germany 


A  computer  software  system  called  ASAT  exists  at  MBB  which  allows  an  automatic 
design  of  minimum  weight  structures.  In  this  paper,  the  application  of  this  system  to 
several  structures  is  described. 

It  is  shown  that  a  structural  optimization  system  can  be  very  useful  in  the  prelimi¬ 
nary  design  of  an  airplane,  especially  when  it  consists  of  several  modules  such  as  static 
load  calculation,  deformations  and  stress  calculation  by  finite  elements,  static  aero- 
elastics,  weight  calculation,  unsteady  aerodynamic  forces,  vibration  calculation,  flutter 
calculation,  flutter  and  strength  optimization  which  all  can  be  used  separately  and  in¬ 
dependently  . 

INTRODUCTION 

For  structural  design  of  modern  airplanes,  the  use  of  an  optimization  computer  pro¬ 
gram  is  mandatory  in  order  to  achieve  a  minimum  weight  structure  whilst  taking  into 
account  both  strength  and  aeroelastic  requirements. 

During  a  cooperation  program  (VJ  with  the  U.S.  Air  Force  Flight  Dynamics  Laboratory, 
the  MBB  company  exchanged  several  computer  programs  in  return  for  receiving  the  FASTOP- 
computer-system  /JJ .  This  exchange  took  place  in  1977, and  for  the  last  three  years,  the 
structural  dynamic  group  of  MBB  has  further  refined  the  program  and  also  added  a  static 
aeroelastic  part  to  it.  This  new  system  is  now  called  ASAT.  (Automatische  Struktur-Aus- 
legung  fur  TragflSchen) .  This  paper  deals  with  the  application  of  ASAT. 

Several  structural  examples  are  treated  in  this  paper: 

.  A  simplified  structure  to  show  the  capabilities  of  the  system  (the  analysis  of  this 
structure  was  partly  sponsored  by  the  ZTL-Research  Program  of  the  German  Ministry 
of  Defense) . 

.  Aeroelastic  efficiency  calcs  for  fin  and  rudder. 

.  Structural  layout  of  a  carbon  fibre  composite  Delta  wing. 

TECHNICAL  APPROACH 

The  ASAT-program  is  able  to  size  cantilevered  or  free-free  surface  structures  for 
flutter  speed  or  strength  constraints.  It  is  based  on  a  finite  element  method.  Buckling 
of  elements  is  considered.  Also  minimum  skin  gauges  can  be  a  limiting  factor  for  sizing. 
The  aeroelastic  efficiencies  are  calculated  directly  by  using  the  aerodynamic  influence 
coefficients  -  no  iteration  procedure  is  applied.  The  mathematical  approach  can  be  found 
in  HI  and  £4 7. 

SIZING  OF  A  SIMPLIFIED  METAL  WING  STRUCTURE  FOR  STRENGTH  AND  FLUTTER  CONSTRAINTS 

In  order  to  try  out  the  computer  system,  a  simplified  structural  model  was  chosen 
(Fig.  1).  The  thickness  to  chord  ratio  is  constant  5  %.  The  surface  is  cantilevered. 

The  conditions  which  are  sizing  the  skin  thickness  against  buckling  are  presented 
in  Fig.  2.  Two  aerodynamic  load  cases  were  defined: 


Load  case  1:  Ma  =  0.9,  q  =  5.52 
Load  case  2:  Ma  =  1.4,  q  =  8.28 


(X  =  5.5° 
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the  upper  and  lower  skin  are  printed  for  successive  steps  of  the  optimization  procedure. 
Initially,  a  constant  skin  thickness  is  provided.  After  three  steps  of  the  SOP-module 
(Struktur  OPtimierung)  a  fully  stressed  design  is  reached  where  the  last  weight  change 
is  only  O.T~kg.  This  is  plotted  in  Table  1. 


Weight 

Initial  weight  for 

wassmsm 

constant  skin  thickness 

kg 

51.4 

56.7 

55.0 

54.3 

TABLE  1  Weight  for  structural  optimization  procedure  (SOP) 


In  Table  2,  the  iteration  procedure  is  shown  for  selected  structural  elements.  It 
starts  with  the  skin  thicknesses  of  step  3  of  SOP  and  then  it  iterates  between  FOP 
(Flatter  OPtimierung)  to  fulfill  the  required  flutter  speed  with  a  minimum  weight  in¬ 
crease  and  still  keeps  the  fully  stressed  design  by  running  through  SOP. 


TABLE  2  Optimization  Progress  for  selected  structural  elements 


The  elements  most  important  for  flutter  speed  (stiffness  change)  are  underlined.  It 
is  interesting  to  note  that  for  instance  the  upper  skin  is  mostly  designed  by  strength 
requirements  whereas  the  lower  skin  thickness  can  be  used  to  raise  the  flutter  speed 
by  a  stiffness  change.  After  a  constant  flutter  derivative  for  each  important  flutter 
element  is  reached  then  the  process  is  finished.  Graphically,  this  is  shown  in  Fig.  4. 


Fig.  4  RESULTS  OF  REDESIGN  STUDY 


A  flutter  speed  is  calculated  for  the  initial  fully  stressed  design  (FSD)  being 
700  kts.  After  five  iteration  steps,  the  desired  flutter  speed  of  900  kts  is  reached 
with  an  increase  of  less  than  3  %  of  total  weight.  The  loss  of  flutter  speed  from  2'  to 
2  and  3'  to  3  can  be  explained  this  way:  The  program  uses  the  old  vibration  modes  to 
get  from  1  to  2 ' ,  but  these  modes  change  a  little  which  is  reflected  in  point  2.  When 
the  structural  changes  are  smaller  and  smaller  then  the  modes  stay  practically  the 
same  (see  point  A',  A  and  5).  In  Fig.  5,  the  elastic  deformations  before  and  after  opti¬ 
mization  are  shown. 


In  Fig.  6,  the  vibration  mode  shapes  are  depicted.  From  this  picture,  it  can  be 
seen  that  the  mode  shapes  stay  almost  the  same, and  only  the  frequencies  are  changed. 
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The  program  FOP  has  also  the  possibility  to  increase  flutter  speed  by  mass  balan¬ 
cing.  Seven  mass  positions  at  the  outer  wing  to  apply  balance  masses  were  provided  but 
the  flutter  derivatives  were  so  small  that  this  possibility  was  neglected  automatically. 

The  final  results  are  presented  in  Fig.  8  and  9  as  skin  thicknesses  for  the  upper 
and  lower  skin  before  and  after  optimization.  Also  the  stress  ratios  -  which  should  be 
unity  when  fully  stressed  -  and  the  flutter  derivatives  are  shown. 


FIG.  8  REDESIGN  RESULTS  FOR  UPPER  FIG.  9  REDESIGN  RESULTS  FOR  LOWER 

COVER  SKIN  COVER  SKIN 


The  normal  force  flow  for  load  case  2  is  shown  in  Fig.  10  as  a  typical  example  of  the 
strength  calculation. 
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Fig.  10  NORMAL  FORCE  FLOW  IN  SPAR  DIRECTION  FOR  LOAD  CASE  2 


Comparisons  for  stress  calculations  with  different  elements  are  presented  in  Fig.11. 
This  figure  proves  that  with  a  relative  crude  element  and  mesh  system  good  correlation 
with  analyses  using  more  sophisticated  elements  -  such  as  NASTRAN  (triangular  membrane 
with  linearly  varying  stress)  -  can  be  achieved.  This  is  an  important  result  because  the 
cruder  the  idealization  can  be,  the  less  computer  time  is  needed  to  run  the  optimization 
program. 
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Fig.  11  STRESSES  FOR  DIFFERENT  ELEMENT 
TYPES  AND  MESH  SYSTEMS 


Fig.  12  IDEALIZATION  OF  A  CFC -VERTICAL  TAIL 


AEROELASTIC  EFFICIENCY  CALCULATION  FOR  A  CFC  FIN  AND  RUDDER 

For  the  structural  desion  of  fin  and  rudder, stiffness  considerations  are  over¬ 
riding  and  not  strength.  Also  the  size  of  these  surfaces  is  influenced  by  their  aero- 
elastic  efficiency.  The  program  ASAT  was  used  to  calculate  the  efficiency  for  a  fin 
and  rudder  of  a  modern  fighter  plane.  The  aeroelastic  deformations  were  calculated 
directly  without  using  any  iteration  procedure  which  is  possible  because  a  full  matrix 
of  aerodynamic  influence  coefficients  is  produced  by  the  aerodynamic  module  of  ASAT. 

The  properties  of  CFC  were  introduced  by  the  stress-strain  law. 

Fig.  12  shows  the  structural  idealization  for  a  CFC-fin  and  rudder.  Fig.  13  shows 
the  deflections  of  the  structure  due  to  a  steady  load  case. 


SEMI -SPAN  :  1990  MM 


TOTAL  LOAD 


Pu  =  98.5  KN 


,  MAX.  DEFLECTION  : 

\ 

UNDEFORMED  STRUCTURE 


Fig.  13  DEFLECTION  DUE  TO  A  STATIC  LOAD 

Fig.  14  presents  relatively  large  changes  in  the  pressure  distributions  due  to 
elastic  fin  deformations  especially  for  the  subsonic  case. 
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Fig.  14  PRESSURE  DISTRIBUTIONS  AND  DEFOR-  Fig.  15  PRESSURE  DISTRIBUTIONS  AND  DEFOR¬ 
MATIONS  AT  SECTION  A-A  DUE  TO  A  MATIONS  AT  SECTION  A-A  DUE  TO  A 

FIN  ANGLE  OF  ATTACK  /3  RUDDER  ROTATION  ANGLE  6 
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Fig.  IS  is  depicting  even  larger  effects  on  the  fin  pressure  distribution  due  to 
rudder  angle  when  elastic  effects  are  considered. 

Fig.  16  and  Fig.  17  present  the  aeroelastlc  effectiveness  of  fin  and  rudder,  and  it 
is  shown  that  the  requirements  which  were  postulated  by  the  aerodynamic  department  can 
almost  be  fulfilled. 


Fig.  16  VERTICAL  TAIL  AEROELASTIC  Fig.  17  RODDER  EFFECTIVENESS 

EFFECTIVENESS 


STRUCTURAL  OPTIMIZATION  OF  A  CFC-DELTA-WING 

For  a  preliminary  design  of  a  CFC-Delta-Wing ,  a  structural  optimization  was  perfor¬ 
med  to  achieve  a  minimum  weight  structure  by  retaining  sufficient  control  surface 
effectiveness.  An  additional  constraint  -  a  certain  amount  of  wing  twist  off  at  a  high 
g-manouvre  -  had  to  be  fulfilled  as  well  £|7. 

The  direction  of  the  laminates  were  selected  in  prestudies  by  the  MBB-stress  de¬ 
partment  which  have  accumulated  a  lot  of  experience  with  the  CFC  material  over  the  last 
years  ^§7. 

Despite  the  fact  that  a  huge  amount  of  papers  has  been  published  lately  about  aero- 
elastic  tailoring  with  CFC,  it  is  our  opinion  that  the  possibilities  for  laying  the 
laminates  are  limited  for  two  major  reasons: 

.  Material  properties  are  only  known  for  specific  composites. 

.  Production  considerations  are  dominating. 

For  these  reasons,  we  took  the  preselected  material  properties  (1J  and  fed  it  into 
the  ASAT  program  as 


Ca  x 

! 

An 

A1 2 

A1  3 

£x 

d>  y 

_ 

A21 

A22 

A23 

r  xy. 

A31 

A32 

A33  . 

V  xy 

DETERMINATION  OF  LOADS 


Two  load  cases  were  chosen  according  to  the  definitions  of  our  loads  group  ^§7: 

Load  case  1 :  Symmetrical  high  g  pullout 

This  is  a  manouvre  case  in  the  subsonic  regime 
giving  the  highest  shear  force  and  bending  moment 
at  the  wing  root. 

Load  case  2:  This  is  a  roll  case  in  the  supersonic  flight  regime 

were  the  aerodynamic  center  of  pressure  is  back¬ 
ward.  This  case  is  not  symmetrical  (initiated  by 
the  ailerons)  but  was  applied  to  both  wings 
symmetrically. 

In  order  to  make  the  loads  calculated  by  ASAT  comparable  to  the  loads  from  our  loads 
group,  the  wing  attitude  and  aileron  angles  are  somewhat  different  -  the  presence  of  a 
canard  had  to  be  reflected  which  the  ASAT  program  is  not  able  to  consider  at  the  moment. 


RESULTS  OF  DEFORMATION  CALCULATIONS  WITH  ASAT  AND  COMPARISONS 

After  establishing  the  structural  model  and  the  loads,  deformation  calculations  were 
performed  which  match  the  stress  group  results  very  well.  Implicitly,  this  is  also  a 
prove  that  comparable  loads  were  used.  Fig.  21  shows  the  vertical  deflection  along  the 
wing  span  for  load  case  1 . 


SrAiVMSE  DIRECTION  Y  [m] 


Fig.  21  COMPARISON  OF  WING  DEFLECTIONS 


In  Fig.  22,  the  wing  twist  angle  along  the  span  is  depicted.  The  4°  twist  off 
angle  at  the  wing  tip  fulfills  the  requirement  coming  from  aerodynamic  performance. 


Fig.  22  COMPARISON  OF  WING  TWIST  ANGLE  04 
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In  Fig.  23,  the  internal  skin  load  distribution  for  load  case  1  is  presented.  From 
this  figure,  it  can  be  deduced  that  extreme  care  must  be  taken  to  accomodate  such  high 
local  forces  into  the  two  rear  CFC  attachments  fitted  to  the  fuselage. 


Fig.  23  INTERNAL  SKIN  LOAD  DISTRIBUTION  Fig.  24  WEIGHT  OPTIMIZED  SKIN  THICKNESS 

ACCORDING  TO  STRENGTH  CRITERIA  (SUB 
SONIC  AND  SUPERSONIC  LOAD  CASES) 

RESULTS  OF  THE  OPTIMIZATION  PROCESS 

After  three  steps  of  the  SOP  program,  the  skin  weight  stayed  almost  constant.  The 
final  weight  was 


Step 

Skin  Weight  ^kg7 

1 

221  .7 

2 

164.4 

3 

163.9 

4 

164.1 

practically  reached  after  the  first  step  but  the  convergence  had  to  be  proven. 

The  weight  saving  amounts  to  about  5  %  of  the  total  wing  weight  which  is  a  very 
considerable  achievement. 

Fig.  24  shows  the  skin  thickness  distributions  before  and  after  structural  optimiza¬ 
tion. 
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DYNAMIC  PRESSURE  [kN/m*] 


Fig.  25  WING-FLAP  CONTROL  EFFECTIVENESS 


The  aeroelastic  efficiencies  for  the  aileron  are  shown  in  Fig.  25  together  with  the 
effectiveness  definition  to  fulfill  the  roll  requirement.  For  the  supersonic  case  where 
the  roll  manouvre  is  iniated  with  the  two  inner  flaps  only  we  have  higher  than  required 
efficiency.  For  the  subsonic  pullout  manouvre,  the  efficiency  for  all  four  flaps  is 
somewhat  below  the  requirement  but  it  still  is  sufficient. 


CONCLUSIONS 

In  this  paper,  it  was  shown  that  the  very  useful  structural  optimization  program 
ASAT  exists  at  MBB  which  was  used  for  several  practical  design  studies. 


The  major  advantage  of  the  system  is  that  it  merges  several  airplane  designing 
disciplines  such  as: 

.  static  loads 
.  stress  calculations 
.  unsteady  aerodynamics 
.  flutter  calculations 
.  static  aeroelastics 
.  weights 


For  this  reason,  communication  errors  are  avoided. 


Due  to  the  versatility  of  the  computer  system,  separate  modules  of  it  can  be  used 
solely,  and  it  is  also  possible  to  make  cross  checks  with  results  from  other  structural 
design  groups.  CFC  structures  can  be  treated  efficiently,  and  the  design  goals  postula¬ 
ted  from  aerodynamic  performance  could  be  reached. 
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A  computer  software  system  called  ASAT  exists  at  MBB  which  allows  an  automatic  design 
of  minimum  weight  structures.  In  this  paper,  the  application  of  this  system  to  several 
structures  is  described. 


It  is  shown  that  a  structural  optimization  system  can  be  very  useful  in  the  preliminary 
design  of  an  airplane,  especially  when  it  consists  of  several  modules  such  as  static  load 
calculation,  deformations  and  stress  calculation  by  finite  elements,  static  aeroelastics,  weight 
calculation,  unsteady  aerodynamic  forces,  vibration  calculation,  flutter  calculation,  flutter 
and  strength  optimization  which  all  can  be  used  separately  and  independently. 
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